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Abstract 

We study the stability of non-Abelian semi-local vortices based on an N = 2 supersym- 
metric H = SU(N C ) x f/(l)/Zjv c ~ U(N C ) gauge theory with an arbitrary number of flavors 
(Nf > N c ) in the fundamental representation, when certain M = 1 mass terms are present, 
making the vortex solutions no longer BPS-saturated. Local (ANO-like) vortices are found 
to be stable against fluctuations in the transverse directions. Strong evidence is found that 
the ANO-like vortices are actually the true minima. In other words, the semidocal moduli, 
which are present in the BPS limit, disappear in our non-BPS system, leaving the vortex 
with the orientational moduli CP^" 1 only. We discuss the implications of this fact on the 
system in which the U(N C ) model arises as the low-energy approximation of an underlying 
e.g. G = SU(N C + 1) gauge theory. 



1 Introduction 



The last several years have witnessed a remarkable progress in our understanding of physics 
of solitons, in particular, that of vortices in non-Abelian gauge theories. These problems are 
important as they are intimately related to one of the deepest issues of particle physics such 
as confinement pQ; they may be important in some condensed- matter physics or cosmological 
problem as well. The progress has been made mainly in the context of supersymmetric gauge 
theories where many aspects of calculations such as the quantum modification of the potential and 
symmetry realizations, are under much better control. After the discovery of vortices carrying 
continuous, non-Abelian moduli in [21 E], many related studies have been carried out [U El El 
El El El HD] • For good reviews, see Refs. [TH [T21 EH E] . Some recent developments involve the 
following issues: 

• A detailed study of moduli and the transformation properties of higher winding vortices 
[T5l fT6l im [T8] has been performed. The last of these papers contains basic results on the 
vortex transformation properties which allow them to be interpreted in a simple group- 
theoretical language. 

• Another development concerns semi-local vortices in systems with larger number of flavors 
(matter fields) [21 [191 EH]- As it happens in the U(l) Higgs systems with more than one 
charged scalar field [2H E2], a new kind of moduli emerges; the thickness of the flux tube 
is no longer a fixed quantity, but corresponds to one of these moduli. This makes the 
corresponding vortex moduli space very interesting. In particular, a new type of (Seiberg- 
like) duality has been found, between different models having dual vortex moduli (and 
sharing the common sigma- model lump limit at very strong couplings) [20] . 

• A systematic study of non-BPS local vortices has been initiated in Refs. [23], El] where 
several kinds of perturbations have been considered. The interactions have a very rich 
structure, depending both on the distance and relative orientations, a new feature not 
present in the Abelian counterpart. 

• A particularly significant development concerns the extension of the analysis to systems 
with generic gauge groups [23 EH1 EI]- As compared to the U(N) models studied in most 
papers, the systems based on, e.g. U(l) x SO(N), U(l) x USp(2N) are characterized by a 
larger vacuum degeneracy, even after the "color-flavor locked" vacuum has been chosen to 
study the solitons. 

• The question of dynamical Abelianization has recently been addressed more carefully in 
Ref. [28]. Only under definite conditions, the non-Abelian moduli fluctuating along the 
vortex length and in time (which are described by a two-dimensional sigma- model), are 
absorbed in the monopoles, turn into the dual gauge group, before Abelianizing dynami- 
cally. The recent result on the non-Abelian vortices with a product moduli [29] seems to 
be particularly significant in this context. 

• Further important developments have been achieved in the study of BPS non-Abelian 
vortices in theories with M = 1 supersymmetry (heterotic vortices) [301 EH] • These works 
generalize the known links between the vortex and the bulk theory to a less supersymmetric 
case. It is reasonable to think that heterotic vortices play a role in the physics of Seiberg 
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duality [ST] . However, so far, there seems to be no clean statement about this. Some related 
works can be found in Refs. [32] . 

• Another direction involves the non-Abelian vortex in the Higgs vacuum of M = 1* theory 
with gauge group SU(N C ), which is a mass deformation of the conformal M = 4 theory. 
This issue has been first studied for N c = 2 in Ref. [33]. In Ref. [31] the case of larger N c 
has been studied, both in the weakly coupled field theory and in the IIB string dual (the 
Polchinski-Strassler background [35]). 

The present investigation is a natural extension of these lines of research, in particular along 
the second and third. Namely, we study the stability of the semi-local vortices in the presence of 
an M = 1 perturbation which makes the vortices non-BPS. In the case of a theory with a mass 
gap, the 't Hooft standard classification of possible massive phases [36] is applicable; as a result, 
there is a very clear relation between the Higgs mechanism and magnetic confinement. In the 
model discussed in this paper, this is not a priori clear due to the presence of massless degrees 
of freedom (see Ref. [37] for some of the subtleties associated with the massless case). 

The question is important from the point of view that the vortex system (with gauge group 
H) being studied is actually a low-energy approximation of an underlying system with a larger 
gauge group, G. The homotopy-map argument shows that the regular monopoles arising from the 
symmetry breaking G — > H are actually unstable in the full theory, when the low-energy VEVs 
breaking completely the gauge group are taken into account. In other words, such monopoles 
are actually confined by the vortices developed in the low-energy H system: this allows us to 
relate the continuous moduli and group transformation properties of the vortices to those of 
the monopoles appearing at the ends, thus explaining the origin of the dual gauge groups, such 
as the Goddard-Nuyts-Olive- Weinberg (GNOW) duals [3B]- The stability problem on non-BPS 
non-Abelian strings have also been studied in another system: the Seiberg-dual theory of the 
M =1 supersymmetric SO(N c ) QCD [39]. 

In these considerations a subtle but crucial point is that when small terms arising from the 
symmetry breaking G — > H are taken into account, neither the high energy system (describing 
the symmetry breaking G — ► H and regular monopoles) nor the low-energy system (in the Higgs 
phase H — > 1 describing the vortices) is BPS-saturated any longer. This on the one hand 
allows the monopoles and vortices to be related in a one-to-one map (each vortex ends up on 
a monopole); on the other hand, the system is no longer BPS and we must carefully check the 
fate of the moduli of the BPS vortices which do not survive the perturbation. In fact, zero 
modes related to global symmetries (orientational modes) still survive, because of their origin as 
Goldstone bosons, while other modes are no more protected by supersymmetry. These issues are 
the subject for investigation in this paper. 

2 Stability of Non-Abelian Semi-local Vortices 
2.1 Semi-local Vortices 

The question of stability of the semi-local vortices [21], [22] for Abelian systems has been investi- 
gated by Hindmarsh [30] and other authors [3H [32] , some time ago. The model which has been 
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studied is an Abelian Higgs system with more than one flavor (iVf > 1), 

c = + ^0(^0) t - \ (0 - e) 2 , (2.1) 

where T>^ = — iA^ is the standard covariant derivative, <fi = (0i, 02; • • • , 4>n { ) represents a set 
of complex scalar matter fields of the same charge. This model is sometimes called the semi-local 
model since not all global symmetries, i.e. U(Nf) here, are gauged. As a consequence, the vacuum 
manifold M is CP^" 1 = SU(N t )/[SU(N f - 1) x 17(1)]. Since the first homotopy group of M 
is trivial for TVf > 1, vortex solutions are not necessarily stable. For (3 = A/e 2 < 1 (i.e. type I 
superconductors) the vortex of ANO-type jHHH] is found to be stable. In the interesting special 
(BPS) case, (3 = 1, there is a family of vortex solutions with the same tension, T = 2tt^. Except 
for the special values of the moduli (in the space of solutions), which represent the ANO vortex 
(sometimes called a "local vortex"), the vortex has a power-like tail in the profile function, hence 
the vortex width (thickness of the string) can be of an arbitrary siz^]. In the limit of large size, 
the vortex essentially reduces to the CP Nf ~ l sigma- model lump (or two-dimensional skyrmion), 
characterized by ^(CP /_1 ) = Z. 

For (3 > 1 (i.e. type II superconductors), vortices are found to be unstable against fluctuations 
of the extra fields (flavors) which increase the size (and spread out the flux). 

Properties of non-BPS solutions in the case of non- Abelian vortices, within the models similar 
to Eq. (12.11) but with a U(N C ) gauge group, have been studied by some of us [21]. The interactions 
among the vortices are found to depend on the relative orientations carried by these non- Abelian 
vortices as well as the distances between the vortices. The two new regimesj called type I*/II* 
have been found in addition to the usual type I/II superconductors [21]. 

The same authors investigated furthermore another class of (type I/I*) systems in a super 
Yang-Mills theory as well |23| . characterized by certain non- vanishing adjoint scalar masses. 
These models are potentially important as they are exactly the sort of systems arising as the 
low-energy approximation as a result of symmetry breaking at some higher mass scale. Under 
such circumstances, the properties of the vortices in the low-energy system are closely related to 
those of the regular monopoles arising at high energies. With this motivation in mind, we shall 
here concentrate on this class of non-BPS vortices. 



2.2 The Model and the Vortex Solution 

The model we consider is an M = 2 supersymmetric (U(N C ) gauge theory with TVf > iV c hyper- 
multiplets (quarks) Qf,Qf (/ = 1, . . . , TVf). Since all the essential features are already present 
in the minimal case i.e. N c = 2, we will concentrate in the following on this case, i.e., the U(2) 
mo del The superpotential at hand is 



V2 



Qfiao + diT 1 + V2m f )Q f + W e (a ) + W g (a i T i )] , (2.2 



*This type of vortex solutions has been termed "semi-local vortices". Though it is not an entirely adequate 
term, we shall stick to it, as it is commonly used in the literature. 

^At large distance, the type I* interaction is attractive for parallel vortices and repulsive for anti-parallel 
vortices. Vice versa for type II* interaction. 

* The extension to more general cases, N c > 2, will be discussed in appendix lAl 
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where 



W e = -f a + 77 a 2 , 



(2.3) 



Two real positive mass parameters 77 and /1 have been introduced for the adjoint scalars. The 
terms proportional to 77, /i break M = 2 supersymmetry (SUSY) to M = 1. mj are the (bare) 
quark masses. £ is the F-term Fayet-Iliopoulos (FI) parameter. In fact, it is S'[/(2) j R-equivalent 
to the standard FI term. If only the £ term is kept, the system remains M = 2 supersymmetric. 

This kind of system naturally arises in the M = 2, SU (3) SQCD, with SUSY softly broken 
down to M = 1 with a mass term of the form 



W = kW 



(2.4) 



Indeed, when the bare masses for the squarks are tuned to some special values, there exist 
quantum vacua in which the non-Abelian gauge symmetry £77(2) x U(l) is preserved (the so- 
called r = 2 vacua) [35]. The low energy effective theory in these vacua is exactly the theory 
we are studying here. As the quantum vacua with an SU(2) magnetic gauge group require 
the presence of a sufficient number of flavors, Nf > 2r = 4, we must necessarily deal with a 
system with an excess number of flavors (Nf > N c ). As these systems in the BPS approximation 
(77, fi = 0) contain semi-local vortices with arbitrarily large widths, the necessity of studying the 
fate of these vortices in the presence of the perturbations 77, \x presents itself quite naturally. 

The bosonic part of the Lagrangian is (we use the same symbols for the scalars as for the 
corresponding superfields): 



+ (VvQfYWQf + VnQfiWQrf - V(Q, Q, a, 



a ) 



(2.5) 



where e is the U(l) gauge coupling and g is the 577(2) gauge coupling. The covariant derivatives 
and field strengths, respectively, are defined by 



n 2 2 



V, (Qf,Q* f ) = U 



Q f ,Q 







(2.6) 



The potential V is the sum of the following D and F terms 



V 

Vi 

v 2 
v 3 



Vx 
9 2 



V 2 
2 

? 



U3 + U4 



e t3k a jak + Tr f [QVQ] - Tr { [Qt 1 Q 



^(Tr f [gtg]_Tr f [gg+; 



|-|Tr f [gr l g]+2/ia ) 



e 
~2 



Tr { lQQ]-£ + 2r)a 



/=i 



1 " r I 21^, 



2^ 



(2.7) 
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where TYf denotes a trace over the flavor indices. The squark multiplets are kept massless in the 
remainder of the paper 



m f = . (2.8) 

The theory has a degenerate set of vacua in the Higgs phase where the gauge symmetry 
is completely broken, which is the cotangent bundle over the complex Grassmanian manifold 
-Mvac = Gr NuNc ~ SU(N { )/[SU(N C ) x SU(N { - iV c ) x U(l)] p] with N c = 2. It is important to 
notice that the first homotopy group of the Grassmanian manifold is trivial (for Nt > N c ). Up 
to gauge and flavor rotations, we can choose the following VEV for the scalar fields 

7e/i oo-o 



Q = Q Vivo i o ■■■ oj ' ao = ' a = ' (2 - 9) 

where a = g^t* /2. The vacuum is invariant under the following global color-flavor locked rotations 
(U c e SU(2), U t e ST/ (2) c SU(N f )$ 

Q -> U C QU} , Q U}QU C , a -> U c aUt , -> U C F^ C ■ (2.10) 

For t] 0, the theory has also another classical vacuum in the Coulomb phase 

g = g f = 0, a = ^, a = 0, (2.11) 

which "runs away" to infinity for rj = 0. In what follows, we consider the vacuum (12.91) and the 
non-Abelian vortices therein. 

Next, we will construct a particular solution for the fundamental (i.e. the minimum winding) 
vortex, simply by embedding the well-known solution for iVf = 2 flavors in our model. In the 
next section the stability of this solution under perturbations of the additional flavors will be 
studied. Setting Q = Q\ the Euler-Lagrange equations for the theory are 

V^F? U = -ig 2 (Q)r{D v Q f - (V^Q^r.Qf) - e ijk (a^a^ + a,Va k ) , 
V»V,Q = —^r, <9^a = -e 2 ^, V »V,a t = -g 2 ^. (2.12) 

We start from the standard Ansatz for the local vortex embedded in the model with additional 
flavors 

n ( Mr)e ie ... \ ,r 3 
Q= { fc(r) ... ) ' fl o = A (r), a = X 1 (r)- , 

M = A? £ = ^ [1 - /i(r)] y , ^ = [1 - /o(r)] , (2.13) 

2 r z 2 r z 

with r 2 = x\ + Notice that the adjoint fields ao, a = ajT*/2 are non-trivial when we consider 
the non-BPS corrections, whereas they vanish everywhere if 77, \x are set to zero (i.e. the BPS 
limit). 



The vacuum is also invariant under pure SU(N{ — 2) flavor rotations which act on the last N{ — 2 columns. 
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By plugging the ansatz ( 12. 13ft into Eq. ( 12. 12ft . we get a set of complicated second order dif- 
ferential equations for the six functions {0 O; 4>ii foy hi A(b Ai} [23J. It is summarized in Appendix 
[A| see Eqs. (1A.14|) - (1A. 19|) with N c = 2. This expression must be solved with the appropriate 
boundary conditions 



A(o) = i, / (o) 

io(oo) = 1, 0i(oo) 



1, /i(oo) = 0, / (oo) = 0, 
1, A (oo) = 0, Ai(oo)=0. 



We find the following behavior for small r 

o ocO(r), lO cO(l), A ocO(l), X t oc 0(1). 
A numerical solution is shown in Fig. [TJ 



(2.14) 



[2.15) 




0.5 1 1.5 



3.5 4 



Figure 1: Profile functions for the vortex in the radial direction r: fa (unbroken/red), f\ (short-dashed/green), 
Ao (long-dashed/blue), Ai (dash-dottcd/black), </>q (dash-double-dotted/magenta), <fii (long-dash short-dash- 
dotted/grey). The numerical values are taken as £ = 2,e = 2,g = 1,?/ = 0.2, /i = 0.4. For this configuration 
T = 0.964T B ps. 



Further solutions can be obtained by acting with the S77(2) c+ f symmetry (I2.10p on this 
Ansatz. As a consequence it develops some internal (orientational) zero modes. In fact, the 
vortex leaves a Z7 (l) c+ f subgroup of S77(2) c+ f unbroken, i.e. these zero modes parameterize the 
space CP 1 = SU{2)/U(1) = S 2 . 

In the BPS limit r\ = \l — 0, the system has vortices which are solutions to the first-order 
equations [3], see also Eqs. ( 1A.23I) and ( 1A.24|) . with tension 



'BPS 



2vre . (2.16) 



As Ao(x) = Ai(x) = (ao = a = 0) in this BPS configuration, its substitution into the Lagrangian 
(12.51) gives the BPS value, T BPS = 27r£, even if it is no longer a solution to the vortex equations 
(12.121) . for ?7,/i 7^ 0. This means that the non-BPS vortex tension derived from Eq. (I2.12p is 
necessarily less than the BPS value. In other words, the model we are considering always yields 
type I superconductivity [23] . 



6 



Our result below - the stability of the ANO-like vortices - thus generalizes naturally the 
known correlation between the type of superconductor and the kinds of stable vortices, found 
in the Abelian Higgs models [101 HTJ H2]. Also, in so much as we study a softly-broken M = 2 
supersymmetric gauge theory and its low-energy vortex system, the present work can be regarded 
as an extension of the one in SU(2) gauge theory [17], even though in the latter the low-energy 
system was naturally Abelian (U (1)). We have verified that the same qualitative conclusion holds 
in that case. 



2.3 Fluctuation Analysis 

In this section, we will generalize the fluctuation analysis of the stability of the Abelian vortex 
of Ref. [ID] to the non-Abelian vortex. As we have mentioned, the first homotopy group of the 
vacuum manifold is trivial, thus the local vortex found in the previous section can be unstable. 
To study the stability of non-BPS local vortices, we must consider the quadratic variations of the 
Lagrangian due to small perturbations of the background fields. It turns out that the quadratic 
variation of the Lagrangian can be written as the sum of two pieces 

5 2 C= 5 2 Cl . ... + S 2 C\ . . .... , (2.17) 

I local fields Iscrm— local fields ' v ' 

where the first term denotes the variation with respect to the fields describing a non-trivial 
background vortex configuration (i.e. the "local" fields), while the second term is the variation 
with respect to the "semi-local" fields (Qf,Qf with N{ > f > N c = 2). The key point is that 
there are no mixed terms (at the second order) between the variations of the background fields 
and the "semi-local" fields. The first term cannot give rise to instabilities (i.e. the local vortices 
with Nf = 2 are topologically stable). Hence, it suffices for our purpose to study only the second 
term 

^Iscmi-iocai fields = WQj'K&Q* + (V h ,5Q 3 ) (V£5Q 3 ) t - , (2.18) 

where we have taken, for simplicity and without loss of generality, N{ = 3. The subscript 'b' 
denotes background fields which we fix to be given by the Ansatz ( 12. 13[) in the following. Let 

us work out this variation explicitly. Even if = Q D in the background fields we keep their 
variations independent. 

Let us first consider the variation of the potential (I2.7p . piece by piece: 



Variation of V\ 

Remember that we are only interested in the variations that involve the semi-local fields, i.e. the 
variation with respect to the third flavor. The crucial observation is that the background value 
for this field is zero: Q^ 3h = = 0. This means that the variation of terms such as TrffQV'Q] 
are already quadratic in the perturbation. Furthermore, the term e^a^ak evaluated on the 
background fields is zero, due to 03 being the only non-zero non-Abelian adjoint field. Thus, the 
variation of V\ with respect to the semi-local fields is at least cubic, coming from the cross terms 
involving the adjoint field. Hence, the quadratic variation vanishes 
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Variation of V2 



Vi contains no adjoint fields. The same argument used for V\ goes through: the variations of V2 
are at least quartic 



5 2 V 2 \ n A = . 



(2.20) 



Variation of V* 



The variation of V 3 gives us the first non-trivial contribution. In this case the variation with 
respect to the semi-local fields is at least quadratic, and it is given by 

,2 



5 2 V 3 \ 
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v Tr f {Q{T 6 Q h ) + 2 /2a 3h ) ( 8Q 3 r 6 5Q 3 + c.c 
+ y (Trt(QlQ h ) - £ + 2 r/aob ) ( WQ 3 + c.c 



(2.21) 



Variation of V4 



The variation of V4 is also quadratic and is simply 



S 2 Va 



Qa,Qs 



\sQ\ (a ob + r 3 a 3b ) 2 5Q 3 + ^5Q 3 (a b + r 3 a 3h ) 2 5Q\ 



(2.22) 



The variations are not diagonal, thus involve mixed terms like 5Q 3 5Q 3 . We can easily diago- 
nalize them by the following change of coordinates (keeping the kinetic terms canonical) 



8Q3 = ~^{q + q ] ) 



5Q 3 = —^-q) 



(2.23) 



which yields the following variation of the potential to second order 

2 

5V L,Q 3 = \ (^f(Qir 3 Qb) + 2/ia 3b ) {\ qi \ 2 - \q T \ 2 - \q n \ 2 + \qii?) 

2 

+ e - [Tr f (Q h Q b ) - £+ 2r/a 0b ) (|g,| 2 - |g,| 2 + \q n \ 2 ~ \hi\ 2 ) 

+ 2 ( a ob + a 3 b) 2 (|<?/| 2 + |g/| 2 ) + g ( a ob - a3b) 2 (|g//| 2 + |<m| 2 ) , 

where the capital indices I and II label the color components. The problem of stability is now 
reduced to studying four decoupled Schrodinger equations. We expand the fluctuations as 



(2.24) 



1i,n = rij; 

k 



(*0 e ikd 



Ql,II 



(fc) 



(2.25) 



Using these expansions, the quadratic variations of the energy density ( 12. 18ft give rise to the 
following Schrodinger equations 



r dr V (ir 



1± (r- \+V 

r dr \ dr ' 



r(k) 



7(*0 



ii>(*0 7(fc) 
m I,IIYI,II ) 



(2.26) 



0.2 1.8 10 
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(a) (b) (c) 

Figure 2: Potentials as functions of r: (a) V} (0) and vff . (b) V} (0) and V} ( 7 0) . (c) Vj {1) and V} (1) . The results are 
plotted for e = 2, g = 1, f = 2, 77 = 0.2 and /ti = 0.4. 

which must be solved with the following boundary conditions at small r: 

4%(r),$%{r) = r k + 0(r k+1 ). (2.27) 

The effective potentials are 

V tii = 4^2 [/o - 1 ± (/1 - 1) + 2 *f + ^(Ao ± AO 2 (2.28) 
2 2 
+ y (0o + <Pi ~ £ + 2 77A0) ± |- (^ - 2 + 2 //Ax) 

Vffi = ^ [/o - 1 ± (A - 1) + 2*] 2 + ^(A ± A1) 2 (2.29) 
2 2 

- y(0o + 0? - f + 2r/A ) T y(0o - 0? + 2//Ai) . (2.30) 

The upper signs refer to the color-index /, while the lower to index //. Since only the first terms 
are dependent on k, it is easy to check that 




-0.05 0.05 0.1 0.15 0.2 0.25 0.3 5 10 15 20 25 30 35 40 

M r 



Figure 3: Left panel: the value of the wave function ip) ' at large distance (r = 20, in the unit of the vortex 
size), as a function of M (see Eqs. (|2.26p ). There are no zeros at negative values of M: this implies the absence 
of normalizable tachyonic modes. Right panel: the wave function for the zero energy fluctuation: M = 0. The 
linear behavior is a common feature for zero energy wave functions. The numerical values of the parameters are 
chosen as in Fig. O 
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(a) 




(b) 




10 15 20 25 30 35 40 



(c) 



Figure 4: Numerical results in the BPS case: e = 2, g = 1, £ = 2, 77 = /a = 0. (a) The potential has a zero 



energy fluctuation, (b) The value of tp\ at a large distance (Y = 20) as function of M. (c) The zero energy wave 
function as function of r. 



for all k. To prove the above inequalities, it is sufficient to recall that the profile functions for 
the gauge fields satisfy: < /o,i < 1. 

From a well-known theorem in one-dimensional quantum mechanics, it follows that (valid for 
the lowest eigenvalues) 



Mf\MP ) 



Mf\M? 



(2.32) 

Hence, in order to exclude the existence of negative eigenvalues, it suffices to study the six 
potentials: v}? and V}? . 

Fig. [2] shows that only the potentials Vf and vjj can be sufficiently negative, such that one 
could expect negative eigenvalues. The potentials Vj and Vj are even divergent at the core 
of the vortex. Indeed, the results of the numerical analysis show (up to the numerical precision) 
that there are no negative eigenvalues for any of the potentials, V®' , V T ' , Vjj and Vfj. We 
have checked this statement in a very wide range of values of the couplings of the theory. This 
strongly suggests the absence of negative eigenvalues for all values of the couplings. The result 
is shown in Fig. [3] for a particular choice of parameters. Note that all the potentials go to zero 
at large distance. This means that there will be a set of positive eigenvalues that represent the 
continuum states, which can be interpreted as interaction modes with the massless Goldstone 
bosons of the bulk theory. 

The status of the zero energy wave functions however is subtler. The wave functions for the 
non-BPS (Fig. [3]) and the BPS cases (Fig. H]) are both non-normalizable, hence they should both 
be interpreted as a part of the continuum. However, there is an important difference between the 
two cases. In the BPS case, the wave function is limited and goes to zero at infinity; its squared 
norm diverges only logarithmically with the transverse volume (~ logL 2 ). In the non-BPS case, 
the wave function is unlimited, and it diverges linearly (Fig. [3]). This implies a divergence of 
its squared norm being quadratic in the transverse volume (~ L 4 ). Such a fluctuation changes 
the vacuum expectation values of the scalars at infinity, thus does not correspond to a size zero 
mode (collective coordinate). The situation is clearer when our vortex system is put into a finite 
volume. In the BPS case, the wave function of the semi-local mode asymptotically approaches 
zero and is an (approximately) acceptable wave function. In a box, it represents a zero energy 
bound state. On the contrary, the wave function for the non-BPS case is non-normalizable. The 
corresponding zero energy state is eliminated and the semi-local excitations have a mass gap. 
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This observation strongly suggests that in the non-BPS case the size moduli disappear. 



As a further check, we note that the potential which gives rise to the size zero modes in the 
BPS limit (which we know exist in that case) is precisely Vj°\ The zero mode is thus given by a 
fluctuation of the field i)f\ Using Eq. (|2~23l) we see that 5Qzi = 5Qsi = v / 2^j°' ) , while all other 
components of Q3 and Q3 vanish. This is exactly the fluctuation we expect for the semi-local 
vortex, which in the BPS limit is known to have the following form [191 120] : 

/ MrY 9 X (r) \ (9 ^ 

W ~ 1 Mr) J ■ {l - 66) 



3 Meta-stability versus Absolute Stability 

The results of the previous section indicate that the local (ANO-like) vortex embedded in a 
model with additional flavors is a local minimum of the action. The problem whether or not this 
solution is a global minimum cannot be addressed with a fluctuation analysis only. In principle, 
the existence of an instability related to large fluctuations of the fields is still not excluded. 

However, we have a reason to believe that the local vortex is truly stable, being the global 
minimum of the action. Note that a vortex cannot dilute in the whole space. In fact, such diluted 
configuration would have an energy that is equal to the BPS bound, = 27r£ = Tbps, while the 
tension of the local vortices, in the class of theories we are considering, is always less: T < Tbps, 
as already noted in Section [2721 To see this, we will first show that the semi-local vortex solution, 
obtained in the BPS limit, is an approximate solution to the non-BPS equations, in the limit 
of an infinite size of the vortex. Then, it is easy to check that the tension of such approximate 
configuration converges to the BPS value. Let us write the equations of motion for the semi-local 
configuration (I2.33P 

e 2 (fM ~<Pl + X 2 ) + Ml + 0i + X 2 ) - 2X 2 ) , 

g 2 (fM + <Pi + x 2 ) + fM - 0? + x 2 ) - 2 X 2 ) , 

^((Ao + AOWA + ^S) , 
^((Xo-X^ + e'A-g'B) , 
^((\ + X i y + e 2 A + g 2 B) , 
e 2 ((Ac + Ax)(0 2 + x 2 ) + (Ac - A x )0 2 + e 2 V A) , 
^ 2 ((Ao + A 1 )(0 2 + x 2 )-(Ao-A 1 )0? + ^5) , (3.1) 
where we have defined 

A = 02 + 02 + x 2_ e + 2r/Ao ? B = (f) 2_ (f) 2 + x 2 + 2fxXi (32) 

For simplicity, we take the gauge couplings to be equal: e = g. Let us consider the following 
set of approximate solutions corresponding to a semi-local BPS vortex in the limit of large size 



X 



f ii _ fo 
/o r 
fi 

0o (fo + /i) 2 0o 



n 



, 0i (fo - hf 

X ' (h + f\-2) 



'X 



K 



Ao 
r 

K 

r 
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(lump limit): a ^> l/ey^: 



(00, 01, X) 




r 2 + \a\ 



fo — fi 



r 2 + \a\ 2 



\ Q = \ 1 = . (3.3) 



The right-hand sides of Eqs. (13 .11) calculated with this configuration vanish obviously. The left- 
hand sides are derivative terms which disappear in the large size limit: the set of fields above, 
trivially satisfy the fourth and the last two equations, while they satisfy the other equations up 
to terms of order (9(l/|a| 2 ) 



fo,i 




8\a\ 2 r 2 


r 




(r 2 + a 2 ) 3 



< 



32 



x" + - 

r 



00 0o (fo + 


hf 




2\a\ 2 r 






V2(|a 


2 _|_ r 2)5/2 


x("2 + /o + 


hf 


V^2(|a 


2|a| 3 


4 r 2 




2 _|_ r 2)5/2 



27 a 



< 



2 (5 

2|a| 2 



Thus, we see that the configuration of Eq. ( 13. 3}) . in the limit \a\ — > 00, is a solution to the 
non-BPS equations of motion: it represents an infinitely diluted vortex, having the tension equal 
to the BPS value. 

We thus see that an infinitely wide vortex is not favored energetically. The only remaining 
possibility is the existence of a true minimum for a vortex with some intermediate size, which 
involves also a non-trivial configuration for the semi-local fields. This case is very unlikely. In 
fact, we could not find any such solution in our numerical surveys. The relaxation method should 
have detected such an unexpected stable vortex. 

We consider our results as a strong evidence for the fact that the local vortex is a true 
minimum, absolutely stable against becoming a semi-local vortex. 



4 Low-energy Effective Theory 

In this short section we only make some considerations on the low-energy effective theory of the 
vortex which follow from the results of the previous sections. 

In the semi-local case, the vacuum has no mass gap, thus there are massless Nambu-Goldstone 
bosons in the bulk. This fact gives rise to subtleties (which we encountered in Sec. 12. 3p when 
defining an effective theory for the vortex only. This is because both the vortex excitations 
and the massless particles in the bulk appear as light modes. As previously mentioned, semi- 
local excitations are non-normalizable and the corresponding states must be considered as bulk 
excitations. In a rigorous effective theory, we should take into account only localizable and 
normalizable modes, which discards the bulk excitations^. This holds both in the BPS and 
non-BPS cases. In fact, effective actions which include semi-local excitations have already been 
proposed in the literature. The effective actions for BPS semi-local vortices derived in Ref. [191 
[TTj . however, are valid only in a finite volume. 

^Similar issues have been already discussed for monopoles and vortices, see for example Refs. [TH H51 
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It seems to us that the following point has not been stressed clearly enough in the literature. 
We consider all the excitations related to non-normalizable zero energy wave functions as part of 
the massless states of the bulk theory. In the BPS case, the surviving supersymmetry enables us to 
describe these fluctuations as collective coordinates of the vortex. It is known, for example, that 
a single non-Abelian BPS vortex has non-normalizable size parameters (in the case Nf > N c ). 
From the point of view that we propose here, however, the effective action on the local vortex 
is given by the orientational part only (normalizable states). The other collective coordinates 
(size) give rise to zero energy fluctuations which represent bulk excitations. 

While these considerations might just be a matter of interpretation in the BPS case, they 
could be important if one wants to push further the study of the non-BPS vortices we started in 
this work. 

5 Conclusion 

In this paper we have examined the stability of the non-Abelian vortices in the context of an 
J\f — 2, H — U(N C ) model with Nf > iV c flavors and an M = 1 perturbation. The particular 
perturbation chosen (the adjoint scalar mass terms) makes our vortices non-BPS. Local (ANO- 
like) vortices are found to be stable; the vortex moduli corresponding to the semi-local vortices 
(present in the BPS case) disappear, leaving intact the orientational moduli - CP Nc ~ l in this 
particular model - related to the exact symmetry of the system. 

This conclusion seems to be very reasonable, as the narrow, ANO-like vortices are needed to 
eliminate the regular monopoles from the spectrum of the full system, in the case our H gauge 
theory arises as a low-energy approximation of an underlying G theory, after the symmetry 
breaking G — > H . If the semi-local vortex would have survived the non-BPS perturbations, the 
magnetic monopoles would be deconfmed [19]. Happily, this is not the case; the Higgs phase of 
the (electric) low-energy theory is actually a magnetic confinement phase. 
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A Extension to U(N C ) 

In this appendix, we will extend the analysis of the vortices in SU(2) x U(l)/Z,2 theory to the 
SU(N C ) x £/(1)/Zjv c with the Nf > N c squarks in the fundamental representation. 

The M = 2 multiplets are the following. The SU(N C ) vector multiplet (A^, a = a± +ia>2), the 
U(l) vector multiplet (A 0fJj , a>o = ao,i + ^0,2) and the hypermultiplets (Q, Q) in the fundamental 
representation of SU(N C ). The squark fields Q and are N c x Nf matrices and their U(l) 
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charges are +1/2 and —1/2, respectively. The covariant derivatives are defined as 
^(Q,Q f ) = (fy - iA„ - L4o M /2)(Q,Q+) , = d M a - z[^,a] . 



(A.1) 



We partially break SUSY by adding mass terms for the adjoint scalars a and a$ in the following. 
The bosonic part of the softly broken M = 2 Lagrangian £ = K — V takes the form 



K 



V 



= Tr 
1 

= Tr 
+ g 2 Ti 



^9 



(A.2) 



.9 



QQ + 2/i(ai + m 2 ) ) ( QQ + 2//(ai + ia 2 



Tr[QQt_Qtg]y 



+ 



Tr [QQ] 



+ 2^(a 0> i + ia , 2 ; 



2Tr 



QQ T + Q'Q )( ai + -a 0>i l N 



. (A.3) 



Our normalization is Tr [T a T b ] = 5 ab /2 for the generators of SU(N C ). The space-time metric is 
taken as rj^ = (+, — , — , — ). The masses of the adjoint scalars (a, a ) are fi and 77. e is the U(l) 
gauge coupling and g is the SU (N c ) gauge coupling. (X) for aniVx^ matrix X stands for the 
traceless part of X0: (X) =X-^^1 N and Tr [(X)] = 0. We will choose the following vacuum 



in which to construct vortex solutions 



Q = Q f 



liv c ,0), 



a = do = . 



(A.4) 



A.l The Vortex Equations 

Let us make an Ansatz for the minimal winding vortex solution. First of all, let us take an 
SU(N C ) generator 



E = diag I 



1 



N c ° V iVc-1 
We make the following diagonal Ansatz 



Nr-1 



Gsu(iV c ), Tr [E 2 



iVc-1 

N r . 



Q = Qi = (diag (e ie Mr), 0i(r), • ■ • , &(r)) , 0) 



^4oi — — Q 



2_^ 

iVr 2 



(1 - / (r)) , a = oq- 



N 



Ao(r) 



^ = -%^(Wi(r))E, a = a 1 = A 1 (r)E 



Note ^QQ + = <QQ f + 2//ai) 



J? + 2/iAi)E holds. 



(A.5) 

(A.6) 
(A.7) 
(A.8) 



II Useful identities: Tr [(X) (Y)} = Tr [XY] - Tr ^ Tr lY] 1 at , (X + Y) = (X) + (Y) 
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Inserting this Ansatz into the Lagrangian above leaves us with an effective Lagrangian for 
the fields /„, fx, A , Ai, O and fa: 



C 

K 

V 



2irr(K - V) , 
JV c -l/( 2 2(AT C -1) V2 2 /, 



N c g 2 r 2 
(JVo - l)9 2 
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N 2 e 2 r 2 N 2 e 2 



N c g 2 

e 2 _4 

JV C ~ ' 2" ' iV c 2 

2 (/ o + (Ar c -l)/ 1 ) 2 0§ + 2(]Vc _ (/ - A) 2 0? 



/ 2 



(A.9) 

2(iV c -l)0 / 1 2 ,(A.10) 



5 2 + ^ 2 + ^ (Ao + (N c - 1)A0 2 + (A - Ai) 2 



JV 2 



_/V 2 r 2 



_/V 2 r 2 



where we have defined 



/vr 4 
2 + (iV c -l)0 2 -^ + ^Ao 



.4 

^ = 0o - 0x + 2/^1 • 
The vortex equations are simply the Euler-Lagrange equations for C: 



./o 



" _ J_o_ 

o 



/i 2^ 



r AT 



[(/o + (N c - l)h)<g + (JV C - l)(/ - fi)<P\] = 
[(/„ + (JVc - - (/o - /i)V 2 ] = , 



r 

a; 



JV c r/e 5 



A + (A + (JV C - l)Ai) 4> 2 + (JV C - 1) (A - Ai) 4>l 



iig 2 B + A (Ao + (JV C - l)Ax) 2 - J- (A - Ai) <t>\ 







+ ^ 
r 



B + -A + — (X + (N c -l)X i y + 



a , (/o + (JV C - 



_/V 2 r 2 



jvT b + J A + jv 2 (Ao " Al) + ^v> 2 ~ 
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. 



These equations should be solved with the boundary conditions 



C/b, fi, A , Ai, 4> , 4>i) 
(/o,/i, A , Ai,0 o ,0i) 



\/f,\/f 



as r — > oo 



(A.ll) 

(A.12) 
(A.13) 

(A.14) 
(A.15) 
(A.16) 
(A.17) 
io = , (A.18) 
(A19) 



(1,1, 0(1), 0(1), 0,0(1)) 



as r — > . 



(A.20) 
(A.21) 



All other solutions can be generated from this Ansatz by a flavor rotation. 

When we turn off the parameters fi and 77, the model recovers full M = 2 SUSY and the 
vortices therein become BPS states. One of the common properties for various BPS states is 
that the EoMs can be reduced to first order differential equations. The energy density can be 
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rewritten as follows 



BPS 



N c -l 



2ixr 



N c g 2 
+ 2 



fx 



- 9 2 W ~ 



I 2 2 


\fo N c e 2 t 


+ N 2 e 2 


r 2 \ 



a fo + (N c -l)h , 

^ iV> 00 



+ 2(iV c - 1) 



^+(iV c -l)0 2 -^ 



/o ~ h 



£— + surface terms 
r 



fo 



(A.22) 



where we have discarded Ao and Ai because of \x = rj = 0. The bound from below is saturated 
for the solutions satisfying the BPS equations 

,2 



/ + (iv c -i)A 

-00 , <Pl 



+ (iv c - 1)0? - ^e) 

fo — h i 



N r r 



N r r 



The tension of the BPS vortex is 



POO 

T= / dr£ BPS = 2n[-f^r)]™ = 2nZ 
Jo 



(A.23) 
(A.24) 

(A.25) 



A. 2 The Perturbations 



In order to study whether the vortex solution in the previous subsection is stable, we perturb the 
fields around the solution considered to be the background. Let us denote the squark fields by 



Q = (Q b + SQ U 8Q 2 ) , Q ] = (Q h + 8Q\,8Q 



(A.26) 



Note 5Qi, SQi are iV c x N c and 5Q 2 , SQ 2 are N c x (N { — N c ). Since the last two are perturbations 
around a vanishing background, they are completely decoupled from all the other fields in the 
quadratic Lagrangian as we have seen in the main body of this paper. To unravel the mixed 
terms, we make the following redefinition of fields 5Q 2 = (q + q^)/y/2 and 5Q 2 = (q^ — q)/y/2 as 
before. In order to derive the usual Schrodinger-type equations, we expand q, q as follows 



ike 



\ 1$L(r) ) 



where ifj^ is an iVf — N c row vector and ip^' is an Nf — N c column vector. 



-ik0 



(A.27) 



Substituting these into the Lagrangian (1A.3I) . we obtain the following four effective La- 
grangians 



£(fc;0) 

2nr 

£(fe;0) 
2-KT 



(dr4%) 2 -Vr(rM K i) 



(*)/ 



,( fc h2 



(*)/ 



\,a) 



jT(k;i) 

27rr 

£{k;i) 

2ixr 



-(dr4 k i) 2 -V^(r)(4 k l)\ (A.28) 
-{dM\? -V«\r){my , (A.29) 
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where A denotes the flavor index: A = 1, 2, • • • , N{ — N c and i = 1, ■ ■ ■ , N c — 1. The Schrodinger 
potentials are defined as 

T ,W_ (^-1)3 2 D| e\ 2(A +(iV c -l)A 1 ) 2 + (jV c - + AT c (fc - 1) >\ 2 



N r 2 N3 V riV, 



1 c 



- (fc) = (N c -l)g 2 B _^ A+ 2(A + (iV c -l)A 1 ) 2 + / 7o + (N c - + iV c (A; - 1) \ 2 ^ 



iV r 2 iV„ 2 V riV, 



c 



As expected, we have obtained only four different Schrodinger-type equations independent of 
N c ,Nf (N c > l,ATf > A/" c ). Therefore, the results we found in the minimal example with N c = 
2, N{ — 3 in the body of this paper are valid more generally Also, when the gauge couplings are 
fine-tuned as e 2 /2 = g 2 /N c , the A^-dependence disappears from the equations altogether. 
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